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SUMMARY 

During dendritic solidification liquid flow is induced both by buoyancy forces and 
solidification shrinkage. Based largely on the experimental base of other NASA- 
sponsored investigators, there is strong evidence that the major reason for the liquid flow 
is the former, i.e., thermosolutal convection. In the microgravity environment, it is 
thought that the thermosolutal convection will be greatly diminished so that convection 
will be confined mainly to the flow of interdendritic liquid required to satisfy the 
solidification shrinkage. 

The major motivation for this work has been to provide improved models to the 
NASA community involved in performing and planning space experiments on dendritic 
solidification with emphasis on convection and macrosegregation. Another motivation 
for this work is that macrosegregation is an extremely important subject to the 
commercial casting community. We believe that work of this type, especially when 
combined with the experimental results of other NASA-sponsored investigators, will be 
crucially important in providing a basic understanding of macrosegregation to the casting 
community. 

This report describes the simulation of thermosolutal convection in directionally 
solidified (DS) alloys. A linear stability analysis was used to predict marginal stability 
curves for a system that comprises a mushy zone underlying an all-liquid zone. In the 
unperturbed and nonconvecting state (i.e., the basic state), isotherms and isoconcentrates 
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are planar and horizontal. The mushy zone is realistically treated as a medium with a 
variable volume fraction of liquid that is consistent with the energy and solute 
conservation equations. The perturbed variables include temperature, concentration of 
solute, and both components of velocity in a two-dimensional system. As a model 
system, an alloy of Pb-20 wt pet Sn, solidifying at a velocity of 2 x 10' 3 cm s' 1 , was 
selected. Dimensional numerical calculations were done to define the marginal stability 
curves in terms of the thermal gradient at the dendrite tips, G L , vs the horizontal wave 
number of the perturbed quantities. For a gravitational constant of 0.01 g Q < g ^ 8 o. 

(g„ = 9.8 cm • s' 2 ), the marginal stability curves show no minima; thus, the system is never 
unconditionally stable. Nevertheless, such calculations quantify the effect of reducing the 
gravitational constant on reducing convection and suggest lateral dimensions of the mold 
for the purpose of suppressing convection. Finally, for a gravitational constant of 10 4 g 0 , 
calculations show that the system is stable for the thermal gradients investigated 
(2.5 <G l < 100 K-cm' 1 ). 

The supercritical thermosolutal convection in directionally solidified dendritic 
alloys was also modeled. The model assumes a nonconvective initial state with planar 
and horizontal isotherms and isococentrates that move upward at a constant 
solidification velocity. The initial state is perturbed, nonlinear calculations are 
performed to model convection of the liquid when the system is unstable, and the results 
are compared with the predictions of a linear stability analysis. The mushy zone is 
modeled as a porous medium of variable porosity, consistent with the volume fraction of 
interdendritic liquid that satisfies the conservation equations of energy and solute 
concentrations. Results are presented for systems involving lead-tin alloys (Pb-10 wt pet 
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Sn and Pb-20 wt pet Sn) and show significant differences with results of plane-front 
solidification. The calculations show that convection in the mushy zone is mainly driven 
by convection in the all-liquid region, and convection of the interdendritic liquid is only 
significant in the upper 20 pet of the mushy zone if it is significant at all. The calculated 
results also show that the systems are stable at reduced gravity levels of the order of 
10" 4 g 0 or when the lateral dimensions of the container are small enough, for stable 
temperature gradients between 2.5 < G L < 100 K-cm" 1 at solidification velocities of 2 to 
8 cm • h' 1 . 

The linear stability analyses and the calculations of the nonlinear thermosolutal 
convection can be used to study the dendritic solidification. In addition to those tasks, 
this report briefly describes analysis and experiments on thermal and salt-finger 
convection in systems comprising a fluid layer overlying a porous medium. Finally, the 
aspect of solute redistribution during dendritic solidification was extended in a model to 
predict the formation of microporosity. Lacking data pertaining to the solidification of 
alloys, methods to estimate necessary transport and thermodynamic properties were also 
developed. 
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I. INTRODUCTION 

The convection of interdendritic liquid in the solid plus liquid region of solidifying 
alloys is responsible for many types of macrosegregation in castings and ingots. The 
modeling of such convection and of the resulting macrosegregation has been the subject 
of many studies since about 1970. However, there are still some important aspects of 
convective phenomena that have not yet been incorporated into dendritic solidification 
theory. In this report, we tie together some fundamental aspects of dendritic 
solidification and of convection in order to model macrosegregation. 

The susceptibility of cast alloys to macrosegregation is dealt with on a day-to-day 
basis by foundrymen, ingot-makers and continuous casting engineers alike. For the most 
part, the approach to solving macrosegregation problems is to evaluate a number of full- 
scale production runs to determine acceptable processing parameters for each new alloy. 
While segregation is sometimes reduced by these empirical methods, it usually persists, 
and often is found to reach unacceptable levels leading to expensive scrap. In many 
processes, macrosegregation in castings, ingots or billets is the overriding factor in 
limiting the size and production rate of the cast product. 

In terms of applications, this report addresses the solidification of directionally 
solidified (DS) castings under conditions in which thermosolutal convection is important. 
During the grant period, we did both linear and nonlinear calculations to predict the 
conditions for the onset of convective instability and to predict the transient nature of 
the convection during the solidification of dendritic alloys. 

Such numerical research can be used in support of NASA experiments designed to 
investigate convective phenomena and macrosegregation in dendritic alloys. In addition 
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to providing an analytical model for the dendritic solidification experiments, numerical 
models for solidification can be used to explore the effect of different orientations of the 
gravity vector on convection and macrosegregation in dendritically solidified alloys in the 
low-g environment. 

n. MACROSEGREGATION AND THERMOSOLUTAL CONVECTION 

A. Early Analyses of Macrosegregation 

In the early works, 11 ' 171 emphasis was on the description of solute redistribution, 
including convection of the interdendritic liquid, and macrosegregation; however, 
thermosolutal convection in the liquid was ignored. It was recognized that because of 
constitutional and thermal effects, the density of the interdendritic liquid varies, spatially 
and temporally within the mushy zone, so convection of the interdendritic liquid was 
predicted and calculated. The first model which included gravity-driven convection of 
interdendritic liquid was done by Mehrabian et al. in 1970 P1 That basic approach, with 

V 

physical and numerical improvements, was subsequently used to model and study 
macrosegregation in cylindrical remelted ingots, 17 ' 9,11 ' 13 ’ 151 under thermal conditions in 
which the width of the mushy zone varied 110,141 and for multicomponent alloys. 110,121 But 
in all cases thermosolutal convection was neglected. 

An example of the quantitative results obtained from one of the early studies is 
shown in Figure 1 in which the convection is relatively extensive. Figure la shows the 
positions and shapes of the liquidus and solidus isotherms in a cylindrical ingot. Because 
the isotherms are concave with respect to the all-liquid region, convection is downward 
and toward the center where the streamlines meet and then point upward (Figure lb). 
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The solute-rich interdendritic liquid flows to the center; consequently macrosegregation 
is strongly positive (i.e., enriched in solute) at that location (Figure lc). In Fig. lb the 
streamlines, of course, do not terminate along the surface of the liquidus isotherm in an 
actual ingot. In the analysis used at the time (1977), convection in the all-liquid zone 
was ignored. Pressure on the liquidus surface was simply given as ambient pressure plus 
the metallostatic pressure, and this was used as a boundary condition for the pressure 
and velocity fields within the mushy zone, itself. Darcy’s law was used to represent 
conservation of momentum in the mushy zone. This treatment is sufficient to describe 
macrosegregation when thermosolutal convection is not important and when the 
geometry of the mushy zone is prescribed. 

Figure 1 is an example in which the segregation is in the form of a continuous 
variation in composition from surface-to-center. Another important manifestation of 
flow of interdendritic liquid is the formation of severe localized segregates, which are 
called "channel-segregates" or "freckles". Freckles were shown (by experiment and by 
calculation) to form when the solidification rate is relatively slow and the interdendritic 
liquid flows in the same direction as and faster than the isotherms. This causes a local 
remelting within the S/L zone and the formation of the channel-type segregates. 17 ' 10,15,161 

B. Thermosolutal Convection 

We emphasize that this reports deals with dendritic alloys. However, practically all 
of the published work pertaining to the analysis of thermosolutal convection in solidifying 
metals deals with "plane-front" processes in which the solid-liquid interface is nominally 
flat. Thus is terms of applications, the available literature is pertinent to crystal growing. 
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Nevertheless, as a starting point to this discussion, reference is made to this body of 
literature. 

It is well known that during the plane front solidification of a binary alloy 
convective as well as morphological instabilities can set in due to thermosolutal 
effects. 118,191 The problem of morphological instability (in the absence of convection) was 
first considered by Tillers et al. [20J and later by Mullins and Sekerka. 1211 It is further 
known that the morphologically unstable plane front may later evolve into a dendritic 
network. 122 ' 241 The evolution of the morphologically unstable solidifying front into 
primary dendrites is qualitatively understood. However, the precise development of the 
dendritic network is still being actively investigated. 

In addition to the morphological instability a convective instability has been shown 
to occur during crystal growth processes; the analysis of this mode was first performed by 
Coriell et al. [251 Their results for Pb-Sn alloys indicate that at low growth velocities the 
preferred instability is, in fact, the convective mode, which could give rise to a nonplanar 
interface and certainly to a solid of nonuniform composition. Their basic analysis was 
recently improved upon, by coupling the convective flow with the morphological 
instability. This was done by Coriell et al. [26] and by Brattkus and Davis. 1271 

The convection which occurs during plane front solidification is the result of the 
interaction of a thermal gradient and a solutal gradient, often deemed doubly diffusive 
convection. 1281 In the literature many studies have been reported on doubly diffusive 
systems showing the development of different types of motions such as finger convection, 
cellular convection, etc. depending upon the directions and the magnitude of thermal and 
the solutal gradients; for a review see Turner. 1291 The research mentioned so far 
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considers the onset of various kind of instabilities. However, the theory is not developed 
well enough to predict the transitions toward the long term behavior of the system. For 
example, it is observed experimentally that during solidification of binary alloys, the 
dendritic network is a common occurrence, and the theories developed for solidification 
with a nominally plane front interface are not applicable to predicting convective 
instabilities in dendritic alloys. 

It is well understood that microsegregation of solute is coincidental with the 
formation of dendrites. These effects are further coupled with convection of the 
interdendritic liquid, and this results in macrosegregation. Macrosegregation can be in 
the form of a freckle or channel-segregate and can occur in a directionally solidified 
casting [30]. Experimental investigations on Pb-Sn and NH 4 C1-H 2 0 systems have 
revealed that the formation of the channel-segregates is sometimes intimately related to 
convection in the all-liquid zones.' 31,32 ' 

Hills et al.' 33 ' developed a full set of thermodynamic equations for the mushy zone. 
Huppert and Worster' 34 ' also formulated a simple mathematical model for the mushy 
zone based on the principles of global conservation relationships. However, their model 
assumes a spatially uniform and temporally constant volume fraction of solid in the 
mushy zone and is not realistic for the case of alloy solidification. In a later study, 
Worster 1351 proposed a model for the development of the mushy zone with time. 

However, his model neglects convection completely, and considers only the mass and 
energy equations. Other studies in the literature on the dynamic behavior of the mushy 
zone include papers by Ridder et al.' 9,13 ' and Maples and Poirier.1 14 ’ 21 ' The primary claim 
of these analyses was to calculate the macrosegregation. Yeum and Poirier' ' analyzed 
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the problem of macrosegregation in a vertical directionally solidified ingot and showed 
that even subtle temperature gradients in the radial direction can cause considerable 
macrosegregation in Pb-Sn dendritic alloys at 1 g 0 (g 0 = 9.8 m • s' 1 ). 

In none of the above works pertaining to dendritic solidification has the problem of 
understanding the onset of thermosolutal convection in the all-liquid zone and/ or the 
mushy zone been attempted. Early in the grant period, a related problem of a doubly 
diffusive system, wherein a porous layer underlies a pool of liquid/ 37 * was studied. 
However, in that work solidification phenomena (i.e., moving mushy zone, heat of 
solidification., etc.) were not included, and the porosity of the porous layer was 
considered to be constant and uniform. Furthermore, the concentration gradients were 
not realistic because the advancement of the mushy zone was neglected (i.e., m 
solidification). 

Particularly noteworthy, in combining solidification and convective phenomena in a 
unified theory, are the publications of Bennon and Incropera who formulated the 
continuity, momentum, solute and energy equations for solid-liquid mixtures/ 38 * They 
used their model to simulate convection and macrosegregation in solidifying aqueous 
solutions of ammonium chloride/ 39,40 * Their computations were for horizontal 
solidification and numerically revealed, for the first time, the formation of channel type 
"A-segregates." 

III. LINEAR STABILITY ANALYSIS OF THERMOSOLUTAL CONVECTION 

We consider directional solidification against gravity, as shown in Figure 2. In 
Figure 2, S denotes the all-solid zone, L is the all-liquid zone, and L+S denotes the two- 
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phase dendritic region, which is usually called the mushy zone. The mushy zone is 
treated as a porous medium with a variable porosity. For the nonconvecting state, we 
assume that the interfaces between the S and the mushy zone and between the mushy 
zone and L are planar. Gravity acts in the negative z direction. 

Figure 3 illustrates the unperturbed solutal and thermal fields in the nonconvecting 
state that are appropriate for the dendritic solidification of an alloy in a eutectic system. 
Corresponding to Figure 2, the mushy zone extends from the eutectic isotherm (z — 0) to 
the isotherm along the dendrite tips (z = z ( ). With the solidification rate, V, less than 
approximately 0.1 cm-s" 1 , which is certainly appropriate in most solidification processes 
of dendritic alloys, we can ignore the effects of curvature at the dendrite tips on both the 
solutal field and the thermal field in the vicinity of z = z r Thus, all isotherms and 
isoconcentrates in the unperturbed state are horizontal. Along the isotherm at the tip of 
the dendrites, the composition C, must be slightly greater than the far-field composition 
C m because there is diffusion in the liquid. Because the Lewis number a/D » 1, the far- 
field temperature T m is beyond the scale of Figure 3 and is not shown, whereas the far- 
field concentration is achieved at approximately z = z, + D/V, where D is the solute 
diffusivity in the liquid and a is the thermal diffusivity. 

We also make the following assumptions that are common to the forms of the 
continuity, momentum, energy, and solute conservation equations that we employ: 

(1) The mushy zone moves with a constant solidification velocity V (in the z 
direction), and the overall shape of the mushy zone is constant. 

(2) There is no mass diffusion in the solid. 
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(3) Flow is two-dimensional. 

(4) The solid and liquid phases have constant and equal densities, p. 

A. The Conservation Equations 

Because solidification is considered to be steady, we introduce a coordinate 
system that moves with the solidification velocity V. The continuity equation, consistent 
with assumptions (3) and (4), is 

- 0 [ 1 ] 

dX dZ 

where u and w are the components of the superficial velocity in the x and z directions. 

For this work, the solid is stationary, and we have selected the momentum 
equation recommended by Beckermann and Viskanta,^ 41 ' with the "effective viscosity" 
taken to be the reference dynamic viscosity, n 0 , and the density of the liquid to be equal 
to the reference density, p Qt except for the body force for which we use the Boussinesq 
approximation. With these assumptions, the components of the momentum equation are 

(x momentum) 


dU 

r r dU 

1 

du 

dU 
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dx 
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In Eqs. [2] and [3], K is the permeability, assumed to be isotropic, and <p is the 
volume fraction of liquid. The reference state is the liquid at a temperature T 0 and 
solute concentration C 0 . The formulation subsumes both the Darcian and non-Darcian 
(Brinkman) terms to account for the porous medium behavior over the entire range of 
volume fraction of liquid within the mushy zone. It is important to note that in the 
vicinity of the dendrite tips, use of the Darcian term alone would introduce serious error, 
whereas by including the Brinkman term, we permit viscous shear. Note that K is a 
function of 0, [42] and K approaches infinity as <p approaches 1. The Boussinesq 
approximation for the buoyancy term is included by the variable density, p, which 
depends on solute concentration and temperature. Finally, v 0 is the kinematic viscosity 

of the liquid at its reference state (C 0 , T 0 ). 

For the energy equation, many additional assumptions are added to those listed 
above. The added assumptions are the following: 

(1) The local temperature of the solid and liquid are equal ( T s = T L = T ). 

(2) The thermal conduction flux is given by Fourier’s rate law for conduction. 

(3) The thermal conductivity within the mushy zone can be represented as an 
effective conductivity, «. 

(4) The latent heat of fusion, L, throughout the solidification temperature range 
is constant. 

(5) It follows from assumption (4) that dH s = dH L = cdT, where H s and H L are 
the intensive enthalpies of the dendritic solid and interdendritic liquid, 
respectively, and c is the heat capacity of the mushy zone. 
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With all of these assumptions, the energy equation can be written as follows: 
(energy) 


il * « il ♦ ( W - v)*L - ivY- k>* . ill* 

at ax 32 £ at c az 


[ 4 ] 


where a = k/ Po c. In Eq. [4], the usual transient, advection, and conduction terms are 
represented; the second term on the right-hand side accounts for the evolution of the 
latent heat of fusion during solidification. 

For the solute conservation equation, we use the model of Flemings and Nereo,^ 
but with no solidification shrinkage and with constant density. The concentration of 
solute in the interdendritic liquid is C, and it is tied to the local temperature by 
equilibrium thermodynamics (i.e., the liquidus of the equilibrium phase diagram). Added 
to the derivation given by Flemings and Nereo,^ 2 ' we include Fickian diffusion of solute 
in the interdendritic liquid in the direction of the thermal gradient. The final form is 
(solute) 


. dC dC , dC n 

<p — + u — + (w - <pV) — - D 

dt dX dZ 


^ + + a* 

dZ dZ dX dX 


- (1 - k) +(1 - k) V — C [5] 
dt dz 


8<p 


where D is the solute diffusivity in the liquid and k is the equilibrium partition ratio, 
which is assumed to be constant, 

Each of Eqs. [1] through [5] reduces to the usual form for the all-liquid zone as 
, <p -* 1. Thus, we are assured of continuity between the variables in the mushy zone and 
the variables in the all-liquid zone. 
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This set of equations requires boundary conditions. They are: 
atz = 0, 

T - T e \ C - C e \ u - w - 0 [6a,b,c] 

at z = oo (far-field conditions), 

C - C m \ T - T n \ u - w - 0 [7a,b,c] 

C E and T e are concentration of solute and temperature at the eutectic isotherm, 
respectively. 

We now express the system of equations in perturbed variables and retain only 
the linear terms. We, therefore, write 


u - u s (z) + u'(x,z,t ) 

[8] 

w - w s (z ) + w'(x,z,t ) 

[9] 

P - p s (z) +p'(x,Z,t ) 

[10] 

C - C s (z) + C'(x,z,t) 

[HI 

T - T s (z ) + T'(x,z,t) 

[12] 

<p - <p s (z ) + <p'(x,Z,t) 

[13] 


where the subscript s indicates the variable in the non-convecting and steady state (i.e., 
the basic state), and the variables u \ w' ,p' , C', T\ and <p' are the perturbed 
quantities. Ultimately, for the linear stability analysis, we solve for the perturbed 
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quantities. Of course, by representing the independent variables as Eqs. [8] through [13], 
it is necessary to describe the basic state. 


B. The Basic State 

By definition of our basic state, u s = w s = 0, and Eqs. [2] and [3] reduce to 


dx 


0 


and 



The solute equation for the basic state comes from Eq. [5]. It is 


[14] 


[15] 


d2c s + 0 - QZ ZZ. c - 0 [16 

dz 2 D dz <p s dz dz <p s D dz * 

The energy equation for the basic state is derived from Eq. [4]. For the mushy zone, it 
is convenient to write it in* terms of concentration, by making use of m = dTjdC , which 
is the slope of the liquidus on the equilibrium phase diagram. The result (0 < z < z,) is 


d 2 C s y dC s VL d< t> s 

dz 2 a dz mca dz 

The terms C s and <p s also satisfy 

C s ( 0) - c £ 
c sN ■ 


[17] 


[18] 

[19] 
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and 


♦.M - 1 ' 201 

Equations [16] through [20] must be solved numerically to obtain C s {z) and <p s (z) 
within the mushy zone. With C s (z) known, T s (z) follows from the liquidus of the 
equilibrium phase diagram. The value of C t is determined as part of the solution for 
concentration in the all-liquid zone. 

For the all-liquid zone, C s and T s are independent, and well-known analytical 
relationships for each are available. For z t < z < «, Eq. [16] reduces to 

£!S ♦ L £S . o pi) 

dz 2 D dz 

The boundary conditions accompanying Eq. [21] are Eq. [19] and 


C s ( «) - C. 


The solution satisfying Eqs. [19], [21], and [22] is 


C s ~ C m 
C, - c m 


- exp 


-Z(z-z,) 
D * 


In order to evaluate C , , we recognize that steady state requires 


dC : 


dz 


[ 22 ] 


[23] 


(C, - C m )V - -D 


[24] 



Also, (dC s / dz ) - G L /m , where G t is defined as the thermal gradient at the tip of the 

dendrites. Then, 


DG, 

C, - C B + 

1 mV 


[25] 


Thus, C s (z ) can be determined for the entire field of 0 < z < °°. Equation [25] is valid, 
provided curvature effects at the dendrite tips are negligible. 

The temperature in the all-liquid zone for the basic state is readily obtained from 
Eq. [4]. It is 


d 2 T t \r dT < 

i + L L - 0 

dz 2 a dz 

The initial conditions for this second-order equation are 


[26] 


and 


f ^ 



[27] 


r,(z,) - T, [28] 

Condition [27] is used because it is convenient to specify G L as a parameter. The 
solution that satisfies Eqs. [26] through [28] is 


T a -T t + G l “ exp 


V 



[29] 
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Thus, T s (z) can be determined for the entire field of 0 < z < «. 


C. The Perturbed Variables 

In addition to the six perturbed dependent variables given in Eqs. [8] through [13], 
it is convenient to define two other perturbed variables. They are 


and 


4> _ 

f-1 

+ 

V 

K 

w 

s 

K 

< j 


[30] 


where, as before, s denotes the basic state, and the variables with primes are the 
perturbed quantities. In reality, this does not add two more dependent variables to the 
problem, because K = K(<p) and p = p(C,T). 

To obtain differential equations for the perturbed quantities, the definitions given 
by Eqs. [8] through [13] and [30] and [31] are substituted into the differential equations 
for continuity, x and z momentum, energy, and solute conservation. For a linear stability 
analysis, the inertial terms in Eqs. [2] and [3] and products of perturbed quantities are 
neglected. In the usual manner, the basic state equations are subtracted from these 
equations. The result of carrying out these operations gives the set of equations for the 
perturbed quantities. These equations and the details of the derivation can be found in 
Nandapurkar et alJ 43 ' 

The far-field conditions for the perturbed variables are 


T'(co,t) - 0 


[32] 
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C'(«m) - o 


[33] 


w'( co.O - i^(oo ,,) - o [34] 

az 

and the conditions at the eutectic isotherm are 

T'(0,0 - 0 I 35 ] 

C'(0,t ) - 0 t 36 ] 

and 


—( 0,0 - 0 [ 37 ] 

dZ 

Equation [37] is applicable for relatively low volume fractions of liquid (e.g., less than 
0.4) for which the permeability is quite low, so that the perturbed velocity can be 
ignored. 

The factor (<p/K) s can be evaluated for the mushy zone with a columnar dendritic 
structure 14 ^ by 




r > 

1 


n - 0 ,] 

K 

i J 

S 

ad\ 

k. 4 


05 


[38] 


where d x is the primary dendrite arm spacing and a is a constant. 

The evaluation of the perturbed density, p ' , is accomplished by recognizing that 
p = P (C,T), so 


p' 


dp Q, + d p y* 

~dC ~dT 
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where (3 C and p T are the solutal and thermal expansion coefficients, respectively. 
Therefore, 

/>' - ~Po[Pc c ' + Pt t '] [39] 

Further, in the mushy zone, where T and C are linked by the equilibrium phase diagram, 
Hi = T jC' , and fi-j- = 

The perturbation equations were simplified by assuming that <p = <p(z), where 
<p{z) — <P S - This makes the numerical solution more tractable, yet maintaining the 
essential ingredient that the stability of the system can be evaluated by perturbing- 

velocity, concentration, and temperature. 

With the assumption that <p ' = 0, it is convenient to treat the all-liquid zone and 
the mushy zone individually and couple them by requiring continuity of the variables and 
their derivatives at z = z, (i.e., at the tip of dendrites). The set of equations for the 
perturbed quantities in the all-liquid zone are simply those used for the mushy zone with 
<p s = 1 and [<p/K\ s = 0. Boundary conditions, at z = « and z = 0, for the perturbed 
variables are given by Eqs. [32] through [34] and Eqs. [35] through [37], respectively. 



The perturbed variables are represented in the usual manner: 


w' - fr(z)e a, * iux 

[40] 

C' - C(z)e ff ‘* iox 

[41] 

T - t(z)e aniux 

[42] 


Here, a is the temporal growth rate, and w is the horizontal wave number of the 
perturbation. By substituting Eqs. [40] through [42] into the conservation equations for 
the perturbed variables, there results a set of equations for which eigenvalues of a satisfy 
the solution. Derivations and the numerical procedures for determining the eigenvalues 
are given in Nandapurkar et alJ 43] 

• To verify the correctness of this methodology, we used it to solve the eigenvalue 
problem studied by Coriell et al. (Eqs. 62 through 66 and 70 through 78 in Reference 25) 
for the case of the directional solidification of an alloy with a planar interface. We 
solved their disturbance equations using our finite difference code, while Coriell et al. 
used a shooting method to obtain the neutral stability curves. In addition, we solved 
their problem by two different approaches. In one, we assumed that the principle of 
exchange of stabilities holds, set a = 0 and computed C, as the eigenvalue. In the other, 
we treated a as the eigenvalue and iterated on C m to obtain the onset of instability 
[. Re(o ) = 0]. Our computed results compared well to the results of Coriell et al. [25] so we 
have confidence in the numerical techniques. 
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D. Marginal Stability Curves 

Marginal stability curves were calculated for Pb-20 wt pet Sn alloy, for various 
fractions of the gravitational constant. Results are shown in Figures 4 through 7. The 
regions identified with S and U represent regions in which the basic nonconvecting states 
are stable and unstable, respectively. The calculations were performed for temperature 
gradients ranging between 0.5 and 100 K-cm \ 

For a gravitational constant of 1 go 0>o = 9.8 ms 2 )> there are t wo branches that 
define the marginal stability of the system (Figure 4). These curves do not exhibit a 
region where the system is unconditionally stable for all wave numbers, so unless the 
horizontal width of the container holding the solidifying alloy is quite narrow, convection 
is expected. However, it would be possible to suppress convection by effecting 
solidification in narrow molds. For example, with G L = 10 K-cm' 1 , the system is stable 
with respect to perturbations with wavelengths less than 0.066 cm. Thus, in molds f 
widths less than approximately 0.66 mm, there would be no convection during 
solidification. 

Figure 5 gives the marginal stability curves for the gravitational constant of 0.5 g 0 . 
Although the nature of the curves here is qualitatively similar to that at 1 go, it is quite 
clear that the region of instability is compressed, as expected. The calculations indicate 
that for this situation, the system would be convectively stable in molds of widths less 
than 0.95 mm for a temperature gradient of 10 K-cm' 1 . 

The marginal stability curves change dramatically with a further reduction in the 
gravitational constant. With the gravitational constant equal to 0.1 g Q (Figure 6), the 
system is not unconditionally stable for any specific value of the thermal gradient. But 
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solidification in wider molds can be accomplished without convection. For example, for 
a gradient of 2 Kcm' 1 , and in molds less than 0.4-cm wide, solidification without 
convection is predicted. Also, for gradients greater than about 2.5 K m' 1 , the system is 
always stable convectively. 

When the gravitational constant is 10' 2 g 0 , the system exhibits two branches 
(Figure 7). The marginal stability branch, drawn with the solid curve, is for stationary 
convection; the region to the left of the broken curve is unstable with respect to 
oscillatory convection. At the reduced level of the gravitational constant, the system is 
stable provided G L > 0.75 K cm' 1 and the wave numbers are greater than approximately 
3.0 cm' 1 . Thus, no convection is predicted for solidification in molds less than 2.0 cm in 
width at 10' 2 g 0 for all gradients greater than 0.75 K cm' 1 . It should be mentioned here 
that when the wave numbers are sufficiently small, the wavelength of the perturbation 
(a = 27 r/w = z) has the same order of magnitude as the dimension of the domain, and 
the assumption of infinitesimality of the perturbations does not hold. 

Finally, the linear stability calculations for a gravitational constant of 10 4 g 0 
showed that the system is stable for all wave numbers, at least for the range of thermal 
gradients analyzed (viz., G L > 2.5 K cm' 1 ). 

E. Discussion 

When analyzing the stability of a solidifying dendritic alloy, the interaction of the 
convection in the all-liquid zone and the convection in the mushy zone must be 
considered. The analysis of the stability of such an interaction, however, appears to be 
quite complex, in that the marginal stability curves exhibit no minima which would define 
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conditions for which the system is unconditionally stable. Further, if an attempt is made 
to describe the stability of the system in terms of the usual thermal and solutal Rayleigh 
numbers, it becomes apparent that these dimensionless parameters are not independent. 
Thus, one cannot impose a constant stabilizing thermal Rayleigh number and study the 
stability of the system as a function of a destabilizing solutal Rayleigh number. 
Unfortunately, this limits the utility of our results to the specific cases considered herein. 
Of course, the same methodology could be applied to other alloys and solidification 
conditions. 

Despite the complex interaction in solidifying dendritic alloys, the conditions for 
marginal stability can be studied in terms of dimensional variables. Such analyses are 
particularly useful in suggesting lateral dimensions of the container that can be used to 
suppress convection at various values of the gravitational constant. Of course, the system 
studied herein was unbounded laterally, so it could be argued that the linear stability 
analysis does not define lateral dimensions to suppress convection. This is true in a 
precise sense; nevertheless, nonlinear calculations of the convecting system confirm the 
predictions of the linear stability analyses.^ Furthermore, as mentioned by McFadden 
et al., [28J the inclusion of the effect of sidewalls in the linear stability analysis for 
thermosolutal convection would make the solution to the problem "intractable, if not 
impossible." However, Guerin et al. l44] state that for the modified problem of predicting 
the onset of convection during directional solidification with a planar interface, the 
presence of the sidewalls is stabilizing. On the basis of their work, we believe that our 
computed results for the onset of convection are conservative (i.e., when our analysis 
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predicts a motionless state, sidewalls would make the system even more stable with 
respect to convection). 

One assumption in our analysis is that the density of the solid equals that of the 
liquid. With unequal densities, there would be a flow toward or away from the solid- 
melt interface. In other words, the basic state would not be motionless, with u s = 0 and 
w s = -cV, where c is the shrinkage. However, because the value of c is small, 
approximately 0.04 for the Pb-Sn alloys analyzed herein, the results would not be 
significantly affected by relaxing this particular assumption. 

Finally, it should be repeated that the volume fraction of liquid, </>, has not been 
perturbed. In a linear analysis, we believe that this is not necessary and is analogous to 
decoupling the morphological stability of the interface from the convective instabilities 
that can occur in solidification with a planar interface. 

IV. ANALYSIS OF NONLINEAR CONVECTION 

* 

In the previous chapter, the thermosolutal convection was analyzed in terms of its 
linear stability. In this chapter, the same type of a system is formulated using 
dimensionless variables, and the supercritical convection is calculated. 

It is important to note that the present nonlinear analysis does not permit the 
volume fraction of liquid in the mushy zone to deviate from its distribution determined 
for the nonconvecting and steady state. While this is a reasonable assumption to make 
for convection near the critical state, and it permits a prediction of the conditions in 
which thermosolutal convection is expected, such an assumption probably does not allow 
an investigation of the supercritical stage of convection that leads to the formation of 
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localized segregates that are sometimes called "channel segregates" or "freckles.''^ 45 ' In 
future work, we intend to model the supercritical convection and the formation of 
channel segregates by relaxing the assumption that the volume fraction of liquid in the 
mushy zone does not vary with time. 

A. Dimensionless Conservation Equations 

The physical situation is illustrated by Figure 2. We assume that a binary alloy 
solidifies vertically in a two-dimensional strip of width W. The coordinate system is 
located at z = 0 and moves upward at the constant solidification velocity, V. 

In the previous chapter, all of the major assumptions are listed. The dimensional 
conservations equations, written for the moving coordinate system are given by Eqs. [1] 
through [5]. We now rewrite Eqs. [1] through [5] with nondimensional quantities. 

(1) The thermal Rayleigh number is 

Ra T - *■ tG _ lH ' 

v 0 a 

where G L is the specified initial thermal gradient at the tip of the dendrites 
and H is a characteristic length. For these nonlinear calculations, we have 
selected H = z„ the height of the mushy zone. 
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(2) The solutal Rayleigh number is 


Rcl c 


gPcC.H 3 


where C„ is the concentration of solute in the bulk liquid. 

(3) The Prandtl number is 

Pr - v 0 /a 

(4) The Schmidt number is 


Sc - v 0 /D 


(5) The Darcy number is 


Da - K/H 2 


(6) The reference velocity is 


U - (gt 3 rG^f* 


(7) The nondimensional interface velocity is 

V - VfU 

The momentum diffusion time scale 
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has been chosen because it lies between the temperature and solute diffusion time scales. 
The temperature and solute concentrations are nondimensionalized according to 


T - 


T' - T 0 
G l H 


and 


C' - C. 

C. 

respectively, where the prime denotes a dimensional quantity and T 0 is the liquidus 
temperature of the alloy with concentration C a . For reference pressure, we choose 

P - P 0 ^ 2 /r 2 

Finally, the nondimensional components of velocity are 


u - u' /U and w - w* /U 


The nondimensional coordinates are 

x-x'/H and z-z'/H 

Nondimensional time is 


t - t'/r 


And nondimensional pressure is 


Notice that in the definition of the nondimensional pressure, a dimensional 
pressure p* is introduced. The dimensional pressure gradient in Eqs. [1] and [2] is 
written as vp' = vp s + vp*, where vp 5 = -p$, so that p* satisfies vp* = vp' + p 0 g. 

Because the equations in the mushy zone reduce to the equations in the liquid 
when 0=1, the whole system can be described by one set of equations. In 
nondimensional form, Eqs. [1] through [5] become 


du dw n 
— + — - 0 

dX dZ 


[43] 


du 

+ 

r > 

Ra T 

1/2 


dU + U du’ 

dt 

Pr 

t 4 


0 

. V J 

dZ 0 dX 


- -0 



dX 



[44] 


dW + 

r i 

Ra T 

1/2 


dW U dW~ 

dt 

Pr 

k. j 


0 

. V J 

dZ 0 dX 


- 0 



— + vV - w 
dz Da 


r ' 

Ra T 

1/2 

RaPrC 
T * c 

Pr 

k A 


Ra T Sc 


[45] 


dT 

dt 



PL d<p 
cG l H dz 


[46] 
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[47] 


3C + 

Ra T 

1/2 

u — + (w - <pp) — 

1 

Uv 2 c + i£ d ± 

dt 

Pr 

K. J 


dX dz 

< J 

Sc 

dz dz 


Ra n 


Pr 


1/2 


f (1 - fc) — (1 + C) 

dz 


B. 77ie Initial State 

We assume that initially there is no convection and that, in the mushy zone, 

Eqs. [16] through [20] of Nandapurkar et al. [43J describe C, T, and <p. The nonconvecting 
state is specified by inputting, G L , V, and C„ with 

<p(x,z t ,0) - 1 

0(x,O,O) > 0 
and 

T'(x, 0,0) - T' e 

In addition to obtaining C, I , and <t> in the nonconvecting mushy zone, the position of the 
dendrite tips, z t , is also obtained as part of the solution. 

The characteristic length, H, is made equal to z t , and the height of the all-liquid 
zone is also made equal to H — z r The temperature at the top of the all-liquid zone 
simply becomes 

T' H -T' t + G L H 
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where T' t is the temperature at the dendrite tips. Thus, within the numerical domain, 
the temperature gradient in the liquid of the initial state is assumed to be linear. Notice 
that this initial state does not exactly match the steady-state thermal field, which is 
strictly obtainable only for an all-liquid zone of infinite height. Because we seek a 
nonlinear calculation for T, C, u, and w, it is not necessary that the initial state must 
match a steady state. However, we choose the values of C, T, and <p that rigorously 
match a steady state for the initial condition of the mushy zone. 

C, Boundary Conditions 

The domain is made finite by means of an artificial top boundary at 
z' = z t + H = 2H and by considering only the region 0 < z < 2 H. 

The boundary condition associated with Eqs. [24] through [27] are as follows: 

(1) No slip condition on velocities along the two vertical walls and the bottom 
boundary. 

(2) No surface tractions at the top open boundary z' = 2H. 

O •/! - 0 

where n is the unit normal vector pointing outward from the boundary 
surface. 

(3) Adiabatic vertical walls. 
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(4) Prescribed temperatures at bottom and top boundaries (with z, - H): 


t(*m - and T(x ' w ' ,) ' * 1 

(5) No solute flux at the vertical and bottom boundaries. 

(6) At the top boundary, the solute must balance the concentration C«, outside, 
viz., at the top boundary, 


_1_ 3C Ra 
Sc dz + . Pr 


VC = 0 


(7) At the bottom, z = 0, 



D. Calculated Results for Nonlinear Convection 

The numerical model uses a finite-element code that is described in Heinrich 
and Heinrich and Yu. [471 The detailed aspects of the method pertaining to a related 
solidification process can be found in Heinrich, [48J and the strategy for selecting a 
nonuniform nodal spacing in the z-direction can be found in Heinrich et alJ 

Calculations were performed for lead-tin systems under a variety of conditions. 
All calculations were started by introducing a random perturbation in the initial solutal 
field at the nodal points of the form 

C(r,z,0) - C 0 (z)(l + 0.005r) 
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where C 0 (z) is the solutal concentration of the initial state and -1 < r < 1, with r taken 
from a random number generator. 

First we discuss the results of nonlinear calculations for the Pb-20 wt pet Sn alloy. 
Table I shows the different cases considered and the results of the linear stability 
analysis. The solidification velocity is 7.2 cm-h" 1 (2 x 10' 3 cm -s' 1 ) in all cases, and the 
upper boundary was chosen at z ' = 2 H. 


Table I. Calculations for Pb-20 Wt Pet Sn 


Calc. 

No. 

Reference 
Thermal 
Gradient 
G l (K cm 1 ) 

8 /go 

(go = 980 cm s' 2 ) 

Container 

Width 

(cm) 

Linear 

Stability 

1 

80 

1 

1.24 

unstable 

2 

80 

10 4 

1.24 

stable 

3 

2.5 

1 

3.46 

unstable 

4 

2.5 

10' 1 

3.46 

unstable 

5 

2.5 

10' 2 

3.46 

unstable 

6 

2.5 

io- 4 

3.46 

stable 

7 

2.5 

1 

0.10 

stable 


In the case of G L = 80 K-cm' 1 , the depth of the dendritic region was found to be 
1.24 cm, and when the gravitational constant isg 0 = 980 cm s' 2 , the non-dimensional 
parameters are Ra r = 4.24 x 10 4 , Ra c = 2.33 x 10 7 , Pr = 1.58 x 10' 2 , Sc = 91.0, and 
Da = 6.15 x 10"V 3 /(1 ' 0)- 

All of our nonlinear calculations agree with the linear stability results in the sense 
that when the stability analysis predicts an unstable case, convection develops in the 
nonlinear system, and when a stable system is predicted, the perturbations die out in the 
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nonlinear model. In Figure 8, the results for calculations 1 after 114 seconds are shown. 

It can be observed that convection is very weak in the all-liquid region, and there is 
essentially no convection in the mushy zone. In fact, the permeabilities used in these 
calculations yield almost no convection in the mushy zone in all cases. 

In the case G L = 2.5 K-cm'\ a depth of 17.3 cm results for the mushy zone, and 
at normal gravity conditions, the nondimensional parameters are Ra T = 4.97 x 10 7 , 

Ra c = 6.28 x 10 10 , Da = 2.75 x 10' V /(1 - <t>) with the same values as before for the 
Prandtl and Schmidt numbers. Results at t = 1080 s for calculations 3 and 5 in Table I 
are shown in Figures 9 and 10, respectively. As expected, much weaker convection is 
observed in Figure 10 and Figure 9 because of the reduced value of g. It can be 
observed that the convection cells in Figure 9 are much longer than those in Figure 8 
because of the weaker stabilizing temperature gradient. These calculations show the 
expected results as predicted by the linear stability analysis of Nandapurkar et alJ 431 
However, these cases show no convection in the mushy zone and hardly any effect of 
thermosolutal convection in the conditions at the dendrite tips. We should point out that 
the linear stability calculations were done for laterally unbounded systems* 4 ^ and for 
values of g sufficiently high to cause convection. The systems were never found to be 
stable unconditionally for all wavelengths. Therefore, we have assumed that a system is 
stable provided that the linear analysis predicts stability for wavelengths equal to or 
greater than the width of the container. 

Another series of calculations was performed for a Pb-10 wt pet Sn alloy at 
G l = 50 Kcm' 1 and a lower solidification velocity of 2.5 cmh' 1 (approximately 
7 x 10 4 cm s' 1 ). These are summarized in Table II. 
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Table II. Calculations for Pb-10 Wt Pet Sn 


Calculation 

Number 

Container Width 
(. H = 2.4 cm) 

8/go 

(g 0 = 980 cm s’ 2 ) 

Linear 

Stability^ 1 

8 

H 

1 

unstable 

9 

2H 

1 

unstable 

10 

H/2 

1 

unstable 

11 

H/ 4 

1 

unstable 

12 

H/ 8 

1 

unstable 

13 

H / 16 

1 

stable 

14 

H 

10" 

stable 


The results for calculations 8 through 12 confirmed that the size of the convection 
cells was independent of the width of the container. The depth of the mushy zone was 
H = 2.4 cm, and the top boundary was chosen at z = 2 H. The parameters for these 
calculations, when the gravitational constant is g 0 , are Ra T = 6.15 x 10 s , Ra c = 9.49 
x 10 7 , Pr = 1.96 x 10' 2 , Sc - 82.3, and Da = 3.947 x 10'V/(1 - <P)- 

Results for calculation 8 in Table II are shown in Figure 11 at / = 700 s. It can 

* 

be observed that, in this case, there is significant convection in the upper part of the 
mushy zone, and convection produces a strong disturbance in the concentration field at 
the tip of the dendrites. The temperature field remains virtually undisturbed. Note that 
dimensionless temperature is defined differently and given in the caption. Figure 12 
shows a magnification close to the dendrite tips at 1887 seconds. A better organized cell 
system can be observed at this later time, and the effect of convection in the upper part 
of the dendrite region is evident. The temperature field has been magnified 100 times 
and shows how weakly it is affected by convection due to the very low Prandtl number. 
In all calculations performed so far, using values generally accepted for the Darcy 
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number, convection in the mushy zone is only significant in the upper 20 pet of the 
dendrite region, if there is any convection at all. Furthermore, as can be observed in 
Figure 12(a), convection in the mushy zone is almost entirely driven by convection in the 
all-liquid region. This is more clearly illustrated by calculation 12 in Table II, for which 
the detail of the flow in the vicinity of the dendrite tips is depicted in Figure 13 at 
1890 seconds. In all these figures, we should keep in mind that the streamline patterns 
obtained are not unique; they depend on the initial perturbations, and, hence, only the 
qualitative features of the flows are significant. 

Further computations have been performed for this system for the parameters of 
calculation 8 in Table II but with the gravitational vector oriented slightly off the 
negative vertical direction. Results show that very small deviations off the vertical 
direction produce a system entirely dominated by natural convection. The effect of both 
thermosolutal and natural convection could only be observed for deviations of up to 
0.01 deg off the vertical direction. For a slightly larger angle of 0.05 deg, natural 
convection completely dominated the flow. These results have been reported 
elsewhere^ 50 ! 

E. Discussion 

Under the simplifying assumptions of steady-state solidification and that the tips 
of the dendrites describe a plane undeformable surface that advances at a constant 
velocity, nonlinear calculations of thermosolutal convection during vertical solidification 
of dendritic alloys have been performed. The model is an extension of the plane-front 
model analyzed by Coriell et alJ 25 ^ and McFadden et alJ 51 ' to include a dendritic zone. 
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The analysis shows significant differences when compared with the plane-front interface 
models, because of the presence of the mushy zone, in that the systems can be stable for 
higher solutal concentration values and were stable for all cases studied at a reduced 
gravity of 10 4 g 0 . This suggests that the plane-front model is not relevant to 
solidification in the presence of a mushy zone. 

The nonlinear calculations are in full agreement with the linear stability analysis 
presented in Nandapurkar et al. l43) and show that the presence of the mushy zone can 
have a stabilizing effect and suppress convection under low gravity conditions or when 
the lateral dimensions of the container are properly chosen. When the system is 
unstable, Tingerlike" convection develops in the all-liquid region that can entrain the 
interdendritic liquid in the mushy zone and significantly affect the solute concentration in 
the neighborhood of the tips of the dendrites. Furthermore, the calculated convection in 
the mushy zone under these conditions is always weak and very sensitive to the choice of 
the permeability function. For the permeability models used in this work, convection in 
the mushy zone never reached further than the upper 20 pet of the mushy zone and 
appeared to be driven mainly by convection in the all-liquid region. 

It was also concluded that at gravity levels of 10 4 g 0 , the Pb-Sn alloys considered 
were always stable to thermosolutal convection; moreover, when gravity is not acting 
exactly vertically, the presence of natural convection fully dominates the dynamics of the 

flow, under the conditions of this model. 

The model presented here is only valid at the onset of convection. After some 
time, constitutional equilibrium is not satisfied, because it is assumed that the 
distribution of the volume fraction liquid remains constant with time. For this reason, 
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calculations have been performed for a short time after the onset of convection that 
describe the main qualitative features of the dynamics of the flow in the vicinity of the 
dendrite tips. Extension of the present model to relax the condition--that the volumetric 
fraction of interdendritic liquid in the mushy zone is constant--is currently being pursued 
with the expectation that it will allow us to model the severe localized segregates in 
which remelting must take place, known as "channel segregates" or "freckles." 

V. OTHER ACCOMPLISHMENTS 

The major emphasis of the research was to model thermosolutal convection 
during solidification of dendritic alloys, and the results of this effort are given in 
Chapters III and IV. In this chapter, we summarize several other tasks of the subject 
research. For the sake of brevity details are not given here, but the many publications 
resulting from the work may be consulted . [37 ' 52 ' 691 

A. Thermal and Salt-Finger Convection in Superposed Fluid and Porous Layers 

When a fluid layer is overlying a porous layer, and fluid in both layers participate 
in the convection process, the nature and characteristics of the convection are quite 
different from those occurring in each layer separately. In directional solidification of 
binary alloys, there is a mushy zone between the solid region and the all-liquid zone. 
Then convection involves both the mushy zone and the all-liquid zone. Of course, in this 
instance, an additional important constraint is the thermodynamics of the phase change 
occurring in the mushy zone. However, in this part of our investigation, we assumed that 


35 


the solid matrix in the mushy zone is inert (i.e., there is no solidification). In this 
manner, rather general results can be obtained based on a small set of parameters. 

We made a linear stability analysis of the onset of thermal and salt-fingering 
convection in the superposed layers. The results show that the presence of the fluid 
layer drastically alters the onset conditions for convection. f37,52J When the thickness of 
the fluid layer is larger than approximately 15% of the porous layer thickness, convection 
is dominated by the fluid layer. When the fluid layer thickness is smaller than 15% of 
the porous layer thickness, convection is dominated by the porous layer. The sizes of 
convection cells in these two cases are drastically different. The prediction of the 
stability theory was confirmed by experiment for the thermal convection case. [53J 

A nonlinear computation was carried out by Dr. Falin Chen, as part of his Ph.D. 
dissertation. [54] The results show that heat transfer is greatly enhanced as the Rayleigh 
number is increased when convection is porous layer dominated. When the convection is 
dominated by the fluid layer, the increase of heat transfer with increasing Rayleigh 
number is rather minimal. Calculations were carried out into the super critical regime 
up to 20 times the critical Rayleigh number. The convection remains steady, and 
oscillations were observed. 

In directional solidification, the solid in the mushy zone is columnar-dendritic, 
so that the mushy zone is a porous medium with anisotropic properties, especially 
in permeability and in heat conductivity. We have recently concluded a linear stability 
analysis of onset of thermal convection in a superposed fluid and anisotropic porous 
layers, and a paper entitled "Convective Instability in Superposed Fluid and 


36 



Anisotropic Porous Layers" by F. Chen and C. F. Chen has been submitted to Physics of 
Fluids. The abstract of the paper is presented below. 

"The onset of thermal convection due to heating from below in a combined 
system consisting of a fluid layer overlying a porous layer whose permeability and 
thermal diffusivity may be anisotropic is studied. The flow in the porous medium is 
assumed to be governed by the Darcy equation, and the Beavers-Joseph condition is 
applied at the interface between the fluid and porous layers. The linearized stability 
equations, together with the boundary conditions, are solved by a hybrid Adams and 
backward difference method. The numerical scheme is first tested for the case in which 
the depth ratio $■ = 10 -4 , where f is the ratio of the thickness of the fluid layer to that of 
the porous layer. The results agree with an exact solution for a pure porous layer to 
three significant figures. Numerical results have been obtained to study the effect of the 
depth ratio f, the degree of anisotropy in permeability f (= the ratio of the permeability 
in the horizontal direction to that in the vertical direction), and the degree of anisotropy 
in thermal diffusivity f T (= the ratio of thermal diffusivity in the horizontal direction to 
that in the vertical direction) on the onset of thermal convection. It is found that 
anisotropy of the porous medium has the most profound effect on the onset of thermal 
convection when the depth ratio is small. At f = 0.1, the critical Rayleigh number 
increases as f is decreased from 1.0 to 0.1 at a given value of f x . This is in accordance 
with the $■ = 0 case. However, when f is decreased beyond the critical value of 0.06, the 
critical Rayleigh number decreases with decreasing f and there is a sudden increase in 
the critical wavenumber a m . This sudden change is the result of switching from 
convection dominated by the porous layer to that dominated by the fluid layer. This 
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result is illustrated by streamline plots of the pertinent case. For depth ratios f > 0.2, 
convection is confined within the fluid layer at onset so that the effect of anisotropy 
becomes less important as the depth ratio is increased." 

B. Experimental Study of Directional Solidification of Aqueous Ammonium Chloride 

Solutions 

Systematic experiments on directional solidification of NH 4 C1-H 2 0 solutions were 
carried out to provide ground-based support for the flight experiment "Casting and 
Solidification Technology" CAST with principal investigators M. H. McCay and T. D. 
McCay of UTSI. A detailed examination of the onset of plume convection in the mushy 
layer was made by varying the bottom temperature of the tank, in which the solution 
solidified. The results are reported in a paper entitled "Experimental Study of 
Directional Solidification of Aqueous Ammonium Chloride Solutions" by C. F. Chen and 
F. Chen, which has been submitted to the Journal of Fluid Mechanics. The abstract for 
this paper is presented below. 

"Directional solidification experiments have been carried out using the analog 
casting system of NH 4 C1-H 2 0 solution by cooling it from below with a constant 
temperature surface ranging from -31.5°C to +11.9°C. The NH 4 CI concentration was 
26% in all solutions, with a liquidus temperature of 15 ° C. It was found that finger 
convection occurred in the fluid region right above the mushy layer in all experiments. 
Plume convection with associated chimneys in the mush occurred in experiments with 
bottom temperatures as high as + 11.0°C. However, when the bottom temperature was 
raised to +11.9°C, no plume convection was observed, although finger convection was 
carrying on as usual. A method has been devised to determine the porosity of the mush 
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by computed tomography. Using the mean value of the porosity across the mush layer 
and the permeability calculated by the Kozeny-Carman relationship, the critical solute 
Rayleigh number across the mush layer for onset of plume convection was estimated to 
be between 200 and 250". 

C. Double-Diffusive Convection under Reduced Gravity 

In a paper to appear, 1551 we reported on the effect of reduced gravity on all forms 
of double-diffusive convection. Using property values of lead-tin alloys, it is concluded 
that diffusive instability, whether due to sideways heating or bottom heating, is not likely 
to occur at reduced gravity levels of 10' 2 to 10 4 g 0 , where g 0 is the gravitational 
acceleration at sea level. Salt-finger convection, however, is likely to be generated at 
much reduced gravity levels. This is because the destabilizing effect of the concentration 
gradient is amplified by the Lewis number (the ration of thermal diffusivity to solute 
diffusivity). This has been discussed in depth in a recent AIAA paper. 1561 There is a 
possibility of interaction between salt-finger convection and convection due to surface 
tension gradients when material processing is being carried out in space. Further studies 
need to be carried out (a) to delineate the parameter range in which such combined 
instability may occur and (b) to determine the nature of such convection motions. 

D. Formation of Microporosity during Dendritic Solidification 

A model for the formation of microporosity was included because it is based 
largely on the description of solute redistribution in the mushy zone that is one of the 
important facets in our models. Porosity is a defect faced by casting engineers, so it was 
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felt that our contribution could be useful to the casting community. The transfer of 
these results was by seminar presentations at an industrial laboratory and a university, 
talks at three conferences, and four publications.* 57-601 

Early in the program Al-4.5 wt pet Cu alloy was chosen as a sample alloy. 
Calculations were done for various concentrations of dissolved hydrogen. By calculating 
the partitioning of both hydrogen and copper between the solid and liquid phases in the 
mushy-zone, the concentration of dissolved hydrogen in the interdendritic liquid during 
solidification could be calculated. In turn, this concentration is used to determine the 
pressure of hydrogen gas in the interdendritic liquid according to the following reaction: 

2 H- H 2 (g ) 

where the underline represents hydrogen in the dissolved state. The equilibrium 
constant K for the above reaction is 

K - P g /C 2 H t 48 l 

* 

where 

p g = pressure of H 2 , and 

Ch = concentration of dissolved H in the interdendritic liquid. 

Also 

K - K(T,C,„) [«] 

where 

T = temperature, and 

C cu = concentration of copper in the interdendritic liquid. 
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C H , C Cu and T throughout the mushy zone are calculated during solidification, so 
that P g can be estimated via Eqs. [2] and [3]. Finally, with a knowledge of the surface 
tension (liquid/vapor interface) of the alloy and space available for potential gas-bubbles 
within the interdendritic liquid, the code can be used to predict the formation of 
interdendritic porosity, and ultimately the volume fraction of interdendritic porosity after 
solidification is complete. 

The procedure, briefly described above, was included in the computer-code and 
used to make a number of calculations for the volume fraction of interdendritic porosity. 
Calculated results agreed favorably with data found in the literature. 157,581 

E. Transport and Thermodynamic Properties 

Another important task was to select and organize data for the various 
quantitative models. In particular, data for the densities of the dendritic solid and 
interdendritic liquid, as well as the densities of the eutectic-solid and eutectic-liquid, have 
been reviewed and carefully extrapolated to the solidification temperature range of 
Pb-Sn alloys 1611 and aluminum-rich Al-Cu alloys. 1621 We also did a careful evaluation of 
the enthalpies of the phases during solidification of Pb-Sn alloys. 1631 For the model which 
predicts porosity, a simplified theory for the surface tension of melts of binary alloys was 
developed, 164-661 which was very useful in extrapolating the surface tension data to the 
solidification range. Finally, because of our treatment of enthalpy during solidification, 
we discovered that the heat of mixing effect in alloys can be important in modelling 
solidification phenomena. In this regard, a paper pertaining to directional dendritic 
solidification 1671 and another on morphological stability 1681 were published. 


41 



In connection with separately supported work on the solidification of binary alloys, 
we estimated the temperature and composition dependence of the viscosity for Al-Cu 
alloys,! 69 ' using a relationship based on Andrade’s viscosity equation. 
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Fig. 1-Interdendritic fluid flow and macrosegregation in an ingot of Sn-15 wt pet Pb. 
(a) liquidus and solidus isotherms; (b) flow lines of interdendritic liquid, 

(c) macrosegregation. From Kou et al.^ 
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Fig. 2--Directional solidification of a dendritic alloy with a mushy zone of depth z t . 



Fig. 3--Solutal and thermal fields for the directional solidification of a dendritic alloy. 
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Fig. 4--Marginal stability curves for the dendritic solidification of Pb-20 wt pet Sn for a 
solidification rate of 2 x 10* 3 cm • s' 1 with a gravitational constant of 1 g 0 . Stationary 
convection is predicted in the region U. 
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Fig. 5-Marginal stability curves for the dendritic solidification of Pb-20 wt pet Sn for a 
solidification rate of 2 x 10' 3 cm ■ s' 1 with a gravitational constant of 0.5 g 0 . Stationary 
convection is predicted in the region U. 
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Fig. 6-Marginal stability curves for the dendritic solidification of Pb-20 wt pet Sn for a 
solidification rate of 2 x 10' 3 cm s' 1 with a gravitational constant of 0.1 g 0 . Stationary 
convection is predicted in the region U. 
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Fig. 7-Marginal stability curves for the dendritic solidification of Pb-20 wt pet Sn for a 
solidification rate of 2 x 10‘ 3 cm s' 1 with a gravitational constant of 0.01 g 0 . The region U 
under the solid curve is for stationary and the region U to the left of the broken curve is 
for oscillatory convection. 
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Fig. 8-Convection in a Pb-20 wt pet Sn alloy (calculation 1) at 114 s with 

G l = 80 Kern' 1 . (a) Isotherms are (T* - T E )/(T, - T E ), (b) solutal isoconcentrates, and 

(c) streamline contours. 
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Fig. 9-Convection in a Pb-20 wt pet Sn alloy (calculation 3) at 1080 s with 

G l = 2.5 Kern' 1 . (a) Isotherms are (T' - T E )/(J t - T E ), (b) solutal isoconcentrates, and 

(c) streamline contours. 
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Fig. 10-Convection in a Pb-20 wt pet Sn alloy of Fig. 9 but at the reduced gravity 0.01 g 0 
(calculation 5). (a) Isotherms are ( T ' - T E )/(T, - T E ), (b) solutal isoconcentrates, and 
(c) streamline contours. 
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16. Abstract 

During dendritic solidification liquid flow is induced both by buoyancy forces and solidification shrinkage. Based largely on the experimental base of other 
NASA -sponsored investigators, there is strong evidence that the major reason for the liquid flow is the former, i.e., thermosolutal convection. In the 
microgravity environment, it is thought that the thermosolutal convection will be greatly diminished so that convection will he confined mainly to the flow of 
interdendritic liquid required to satisfy the solidification shrinkage. The major motivation for this work has been to provide improved models to the NASA 
community involved in performing and planning space experiments on dendritic solidification with emphasis on convection and macrosegregation. Another 
motivation for this work is that macrosegregation is an extremely important subject to the commercial casting community. We believe that work of this type, 
especially when combined with the experimental results of other NAS A-sponsored investigators, will be crucially important in providing a basic understanding 
of macrosegregation to the casting community. This report describes the simulation of thermosolutal convection in directionally solidified (DS) alloys. A 
linear stability analysis was used to predict marginal stability curves for a system that comprises a mushy zone underlying an all-liquid zone. In the unperturbed 
and nonconvecting state (i.e., the basic stale), isotherms and isoconcentrates are planar and horizontal. The mushy zone is realistically treated as a medium with 
a variable volume fraction of liquid that is consistent with the energy and solute conservation equations. The perturbed variables include temperature, 
concentration of solute, and both components of velocity in a two-dimensional system. As a model system, an alloy ofPb-20 wtpctSn, solidifying at a velocity 
of 2 x 10 3 cra-s 1 , was selected. Dimensional numerical calculations were done to define the marginal stability curves In terms of the thermal gradient at the 
dendrite tips, G i , vs the horizontal wave number of the perturbed quantities. For a gravitational constant of O.Olg (g 0 "9.8 cm-s 2 ), the marginal stability 
curves show no minima; thus, the system is never unconditionally stable. Nevertheless, such calculations quantify the effect of reducing the gravitational 
constant on reducing convection and suggest lateral dimensions of the mold for the purpose of suppressing convection. Finally, for a gravitational constant of 
1 calculations show that the system is stable for the thermal gradients investigated (2.5sG^ 100 Kcm 1 ). The supercritical thermosolutal convection in 
directionally solidified dendritic alloys was also modeled. The model assumes a nonconvective initial state with planar and horizontal isotherms and 
isoconcentrates that move upward at a constant solidification velocity. The initial state is perturbed, nonlinear calculations are performed to model convection 
of the liquid when the system is unstable, and the results are compared with the predictions of linear stability analysis. The mushy zone is modeled as a porous 
medium of variable porosity, consistent with the volume fraction of interdendritic liquid that satisfies the conservation equations of energy and solute 
concentrations. Results are presented for systems involving lead-tin alloys (Pb- 10 wtpctSn and Pb-20 wt pet Sn) and show significant differences with results 
of plane-front solidification. The calculations show that convection in the mushy zone is mainly driven by convection in the all-liquid region, and convection of 
the interdendritic liquid is only significant in the upper 20 percent of the mushy zone if it is significant at all. The calculated results also show that the systems 
are stable at reduced gravity levels of the order of 1 0 4 ^ or when the lateral dimensions of the container are small enough, for stable temperature gradients 
between 2.5<G i :£ 100 Kcm 1 at solidification velocities of 2 to 8 cm h A The linear stability analyses and the calculations of the nonlinear thermosolutal 
convection can be used to study the dendritic solidification. In addition to those tasks, this report briefly describes analysis and experiments on thermal and salt- 
finger convection in systems comprising a fluid layer overlying a porous medium. Finally, the aspect of solute redistribution during dendritic solidification was 
extended in a model to predict the formation of microporosity. Lacking data pertaining to the solidification erf alloys, methods to estimate necessary transport 
and thermodynamic properties were also developed. 
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